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ABSTRACT 

We propose an observational test for gravitationally recoiling supermassive black holes 
(BHs) in active galactic nuclei, based on a correlation between the velocities of BHs 
relative to their host galaxies, |At;|, and their obscuring dust column densities, Sdust 
(both measured along the line of sight). We use toy models for the distribution of 
recoil velocities, BH trajectories, and the geometry of obscuring dust tori in galactic 
centres, to simulate 2.5 x 10® random observations of recoiling quasars. BHs with recoil 
velocities comparable to the escape velocity from the galactic centre remain bound to 
the nucleus, and do not fully settle back to the centre of the torus due to dynamical 
friction in a typical quasar lifetime. We find that |Au| and Sdust for these BHs are 
positively correlated. For obscured (Sdust > 0) and for partially obscured (0 < Sdust < 
2.3 g/m^) quasars with |Au| ^ 45 km/s, the sample correlation coefficient between 
logio(|Au|) and Edust is ^45 = 0.28±0.02 and r 45 = 0.13±0.02, respectively. Allowing 
for random ±100km/s errors in |Au| unrelated to the recoil dilutes the correlation 
for the partially obscured quasars to r 45 = 0.026 ± 0.004 measured between |Au| and 
Edust- A random sample of > 3,500 obscured quasars with |Au| ^ 45 km/s would 
allow rejection of the no-correlation hypothesis with Scr significance 95 per cent of 
the time. Finally, we find that the fraction of obscured quasars, Jobs(|Au|), decreases 
with |Au| from J/)bs(< 10 km/s) > 0.8 to J^obs(> 10® km/s) < 0.4. This predicted 
trend can be compared to the observed fraction of type H quasars, and can further 
test combinations of recoil, trajectory, and dust torus models. 

Key words: black hole physics — methods: observational — galaxies: active — 
galaxies: nuclei. 


1 INTRODUCTION 

Both theoretical models fe egelm an. Blandford fc Reed 
Il98d: IV olonteri. Haar dt fc Madaul I 2 OO 3 I ) and observations 
(e.g. IComerford et al.l 20091 ) suggest that it is common for 
supermassive black holes (SMBHs) in galactic centres to 
form binaries that gradually lo se their energy through radi¬ 
ating gravitational waves (e.g. iHaehnefal 19941 : ISesana et al.l 
l2005h . Numerica l simulations of black hole (BH ) binary 
mergers (see e.g. iHealv. Lousto fc Zlochoweill201 j and ref¬ 
erences therein) suggest that the merger remnant SMBHs 
receive a recoil velocity of typically several hundred (and 
in some spin and mass-ratio configurations up to thou¬ 
sands) of km/s, due to highly anisotropic gravitational- 
wave emission in the final merger phase. Taking into ac¬ 
count the spatial distribution of mass in the galactic centre 
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region, and corresponding dynamical friction, this means 
that the SMBH can engage in a damped oscillating mo¬ 
tion llM adau fc Quataert l2004l : iKomossa fc Merritt! 1200^ : 


iTanaka fc Haimanll2009l )with an amplitude comparable to, 
or exceeding the 0(10-100 pc) size of the optically thick 
dusty molecular torus (‘dust torus’) believ ed to be surr ound- 
ing galactic centres (for an overview, see iHonidliool ). For 


a large recoil, the SMBH can escape from the galactic cen- 
tre and can remain wande ring within the dark matter halo 
jGuedes et al.ll2009l. [20TI1 . As accreting material can stay 
gravitationally bound to the moving SMBH, and thus the 
SMBH can remain act ive for 10^ — 10® yr after the merger 
event (e.g. iLoe ^ l2007l) . kinematic and spatial signatures of 
the recoil could be found in the spectra of quasars (QSOs; see 

e^gjBonnmg^Shiekls_^jj^iajr^ij |2007l : iBlecha et al.ll20ril : 


iKomoss^ 2OI2I : Blecha et al.l 20151 ). These signatures could 

confirm the gravitational recoil of merger remnant SMBHs, 
and could probe the SMBH binary parameters, recoil tra- 
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jectory models, as well as the spatial distribution and com¬ 
position of matter in the galactic centre region. 

In this paper, we propose that in the presence of 
dust tori obscuring galactic nuclei, the gravitational re¬ 
coil of active merger remnant SMBHs should introduce a 
correlation between dust column mass densities along the 
line of sight (Edust) and magnitudes of line-of-sight pe¬ 
culiar velocities of SMBHs relative to their host galax¬ 
ies (|Au|). Proxies to estimate both of these quantities 
can be measured from observable fea t ures of QSO spectr a 
dBonning. Shields fc Salviander 2007|;_Ledouj^t a^ 201^ . 


As pointed out earlier bv iKomossa fc Merritt! ( 20081 1 . re- 
coiling BHs can spend a significant fraction of their time 
off-nucleus, possibly reducing the fraction of fully obscured 
(‘type IF) QSOs; we also follow up on this suggestion. 

Using a selected combination of models of gravita¬ 
tional recoil, SMBH trajectories, and obscuring dust tori, 
we demonstrate the feasibility of detecting Edust~|Au| cor¬ 
relations by simulating a set of 2.5 x 10® random observations 
of recoiled QSOs with a Monte Carlo method. We charac¬ 
terize the strength of correlation between Edust and | Au| for 
obscured (i.e. Edust > 0) QSOs in two different |Aw| inter¬ 
vals, and estimate the number of obscured QSOs that could 
be used to reject the hypothesis of no correlation between 
Edust and |Au| with 3 ct signihcance. As the strength of the 
correlation between Edust and |Au|, as well as the underly¬ 
ing Edust (|Au|) relation depends on the presumed combina¬ 
tion of models, observational studies on QSO spectra could 
provide an opportunity for testing chosen combinations of 
these models. We also calculate the fraction of QSOs ob¬ 
scured by their dust tori, J^obs(|Au|), which, compared to 
the observed fraction of type II (i.e. obscured) QSOs could 
provide an independent test for a chosen combination of re¬ 
coil, SMBH trajectory, and dust tori models. We propos e to 
use the SDSS-DRIO Quasar Catalog dParis et al.ll20l3j to 
perform these tests in the near future, and will report on 
this observational search in a separate publication. 

The paper is organized as follows. In § [2] we describe the 
gravitational recoil and SMBH trajectory models we used in 
our Monte Carlo simulations. In § |3l we describe our imple¬ 
mentation of a smooth dust torus model. We present our 
results in § [4] and discuss some issues that will be relevant 
for an observational search for the proposed correlations in 
§[5l Finally, we offer our conclusions and summarize the im¬ 
plications of this work in § [h] 


2 BH DYNAMICS 

Anisotropic emission of gravitational waves emitted by co¬ 
alescing SMBHs carry away linear momentum, resulting in 
a recoil of the merger remnant SMBH in the opposite di¬ 
rection. Numerical simulations of this process have been 
carried out for merging SMBHs with equal and unequal 
masses, zero and non-zero spins aligned or counter-aligned 
with the orbital angular momentum, and with spins point- 
ing in random directions with equal probability (see e.g. 
iHealv. Lousto fc Zlochoweill201 j and references therein). 

To keep our Monte Carlo simulation as general and 
realistic, but at the same time, as computat ionally cheap 
as possible, following iTanaka Sz HaimanI (l2009l ). we adopted 
the analytical formulae given in equation (4) of iBaker et al.l 


1 I 2 OO 8 D to construct the dis tribution of recoil velocity mag¬ 
nitudes, Urecoii (uote that iBaker et al.l I 2 OO 8 I gives results 
very similar to m ore recent, slightly modified formulae in 
iLousto et ahllioi^ and iHealv. Lousto fc Zlochoweij|2014h . 

The direction of the recoil velocity depends on the ori¬ 
entation of the spin and orbital angular momentum vectors. 
It is plausible that prior to merger, the BHs accrete a signifi¬ 
cant amount of gas coherently from a circumbinary accretion 
disc, whose in ner regions are aligned with the binar y’s or¬ 
bital plane fe.g llvanov. Panaloizou fc Polnar"^ JlQQgl lj. One 
then expects that the spin angular momentum vectors at the 
time of the merger may be aligned with the orbital angular 
mom entum of the binary jBogdanovic. Reynolds fc Milleil 
I 2 OO 7 I: although there could still be signihcant misalign¬ 
ment at_jnerge^fOT^at-spinning and unequal-mass bina¬ 
ries; ICerosa et al. IliH). Thus, the kick direction may not 
be random, and may lie preferentially in the plane of the 
circumbinary disc (and also in the binary orbital plane). 
However, we emphasize that here we are interested only in 
the kick direction relative to the orientation of the larger 
scale nuclear torus. It is much less clear whether the bi¬ 
nary orbital plane and/or accretion disc is aligned with the 
symmetry plane of this torus. This depends on the trans¬ 
fer of angular momentum between large (100 pc) and small 
(sub-parsec) scales, which is sens itive to turbulence, star for¬ 
mation, and feedback processes dPubois et al.|[20141 . Obser- 
vationally, parsec-scale maser discs (which could be taken 
as proxies for the accretion discs that determine the kick 
direction) appear to be oriente d randomly with respe ct to 
the plane of their host galaxies llKormendv &: Hal2013ll . We 
therefore simply assume that the kick direction is random, 
i.e. not aligned with the symmetry plane of the torus. Thus, 
we chose directions of recoil velocities randomly from a uni¬ 
form distribution covering a whole sphere. _ _ 

Using the fitting formulae of I Baker et al.l (l2008lj . we 
calculated Urecoii for a given pair of masses (m. 1 , 2 ) and di¬ 
mensionless spin vectors of the merging SMBHs (ai ,2 = 
c5i,2/Gmi 2 , where S'!, 2 are the spins of the SMBHs, c is 
speed of light, and G is Newton’s constant). Observations 
of active galactic nuclei (AGNs) indicate that a signihcant 
number of SMBHs have |qi, 2 | > 0.9, although a second pop¬ 
ulation of SMBHs with 0.4 < |di. 2 l < 0.8 w as found for 
mi ,2 > 4x llF Mq SMBHs 1 see I Revnoldsll 20 1^1 . In our simu¬ 
lations, for simplicity, we hrst assigned each merging SMBH 
to one of the two |5i,2| populations with equal probability, 
and then drew a random |ai, 2 | value from the corresponding 
(i.e. |di, 2 | G [0.9, 1] or |di, 2 | G [0.4, 0.8]) interval. The direc¬ 
tions of both spin vectors were chosen independently from 
a uniform distribution covering a whole sphere. 

We randomized pairs of mas ses independently f r om th e 
SMBH mass function given in lAller fc Richston3 (l2002ll . 
downscaled both mass values by a factor of 2, and used them 
as values for mi and m 2 . The fact that we chose the mass 
of the merger remnant as M = mi -|- m 2 means that in our 
Monte Carlo simulations M values are under-represented at 
the lowest (~ 10® Mq) and at the highest (~ 10® Mq) val¬ 
ues, and overrepresented at value s in between, compared t o 
the SMBH mass function given in lAller fc Richstonel (l2002h . 
Fig.[T]shows the histograms M and UrecoU values we obtained 
by randomizing 2.5 x 10® pairs of merging SMBHs. 

Using the above set of Wrecoii’s with uniformly dis¬ 
tributed random directions, we simulated the resulting 
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Figure 1. Left: histogram of 2.5 x 10^ simulated merger rem¬ 
nant SMBH masses M = mi -\-m 2 (red curves), where the values 
of 2mi and 2m2 were randomi z ed fro m the SMBH mass func¬ 
tion given in lAller & Richstond l|2002l'l and represented by dark 
blue curves. As a result of the randomization process, M val¬ 
ues are under-represented at the lowest 10® Mq) and at the 
highest (~ 10 ®Mq) values, and overrepresented at values in be¬ 
tween, compared to the SMBH mass function. Right: histogram 
of the corresponding recoil velocity magnitudes, Urecoil? of the 
m erger remnan t SMB Hs, obtained using the framework presented 
in lBaker et al.l ll200gt) . and assuming that merging SMBHs have 
dimensionless spin magnitudes, |ai, 2 |, from two disti nct popu¬ 
lation s (i.e. |ai, 2 | ^ [0.9,1] and |ai, 2 | ^ [0.4,0.8]; see iRevnoldsI 
120131) with equal probability, that |oi, 2 | values are uniformly dis¬ 
tributed within each of the two |oi, 2 | intervals, and that ai ,2 
have uniformly distributed random directions. 


SMBH traiectories base d on the model presented in 
iMadau fc QuataertI l|2004h . The initial position of the SMBH 
and the origin of our coordinate system were chosen to co¬ 
incide with the galactic centre. This model assumes that 
the SMBH is embedded in a spherical stellar bulge, and 
is decelerated by dynamical friction. The one-dimensional 
velocity dispersion of stars in the bulge were calculated 
using the empirical M-a relation (|Tremaine et al.l l2002l 'l 
as (tid = (l/v^) X (M /1.3 X 200 km/s . Ac- 

cording to ob s ervat ions bv iBarth. Greene &i Hoi (l2005l ) and 
lOreene fc Hoi (|200fil l. this relation extends to active SMBHs 
with masses as low as M ~ 10® Mq. We simulated the trajec¬ 
tories with a At = 10® yr time resolution up to a randomly 
picked final time between T G [At, 3 x 10“^ At]. The maximum 
duration of Tmax = 3 x 10^ At = 3 x 10^ yr was chosen to 
be compara ble with the m aximum known duration of QSO 
activity (see lMartinill200 J for a review). The simulation was 
terminated before T if the orbit of the SMBH has decayed 
and the SMBH settled back at the galactic centre (i.e. if the 
radial distance, d, and radial velocity, Vr, of the SMBH con¬ 
verged to d < 10“^ pc and Vr < 10~^ km/s, respectively). 
The position, velocity, and Av of the SMBH in its final state 
were calculated and stored, as if they were results of a QSO 
observation made at a random time during the QSO activ¬ 
ity after the recoil. Histograms of d and Vr values obtained 
for the 2.5 x 10® QSOs are shown in Fig. [3 Note that we 
assumed that recoiled SMBHs remain active along their en¬ 
tire trajectory, and also that they are active at the time the 
moc k observation i s mad e. 

lOuedes et al.l (120091 ) followed the motions of merger 
remnant SMBHs with fixed masses in various configurations 
of dark matter haloes using numerical simulations. They 
concluded that due to the asymmetric and inhomogeneous 
mass distribution in the halo, the SMBH trajectories suffer 



Figure 2. Histogram of distance d (in pc; blue curve), and ra¬ 
dial velocity magnitude, Ur (in km/s; red curve), of the 2.5 X 10® 
SMBHs, both measured from the galactic centre at the randomly 
picked time of a mock observation. Both histograms show bi¬ 
modality. The peaks at lower d and Vr values (log,Q(VAL) < — 1) 
correspond to SMBHs that underwent at least a half period of os¬ 
cillation and thus suffered a strong orbital decay due to dynamical 
friction in the galactic centre. The peaks at high (log]^g(VAL) > 2) 
d and Vr values correspond to SMBHs whose initial recoil veloc¬ 
ity was high enough to make them gravitationally unbound and 
escape from the galactic centre region. 


large deviations from the lMadau fc QuataertI l|2004h trajec¬ 
tories. SMBHs might not even return to the galactic centre, 
after the SMBHs reach the first turning point in their os¬ 
cillating motion. To check how much this could affect the 
results of our Monte Carlo simulations, we simulated an in¬ 
dependent set of 10® QSOs, and examined the fraction of 
SMBHs that reach the first turning point before the mock 
observation is made. We found that for the entire QSO sam¬ 
ple, this fraction is ~ 31 per cent, while it is reduced to 
~ 11 and ~ 4 per cent for QSOs with |Au| ^ 5 km/s and 
|Au| ^ 45 km/s, respectively. The reason for the full pop¬ 
ulation having a higher fraction is that once the SMBHs 
make their first U-turn, their orbits decay rapidly due to 
dynamical friction. As a result, many of these end up with 
I An] < 5 km/s and |Au| < 45 km/s at the mock observation 
time (as shown in Fig. [2J. As we will show in § [4] below, only 
QSOs with |Au| ^ 5 km/s show a Edust-jAuj correlation. 
Furthermore, below we propose to study only the subset of 
QSOs with |Au| 1^ 45 km/s, due to the limitations imposed 
by inevitable random velocity errors. We conclude that in 
this sub-sample of QSOs, the effects of asphericity and in¬ 
homogeneities will reduce the predicted correlations only by 
« 4 per cent, and we ignore this complication in the rest of 
this paper. 

A similar deviation from a purely radial trajectory may 
arise on much smaller scales, inside the spherical nuclear 
star cluster (stellar core) enshrouding th e galactic centre. 
As shown bv lOualandris fc MerrittI (|20n8h . when the SMBH 
mass is not negligible compared to the mass of the stel¬ 
lar core, the SMBH and the stellar core exhibit oscillations 
about their common centre of mass. As a result, the SMBH 
is not exposed to the largest dynamical friction near the ge¬ 
ometric centre of the core, and it may take up to ~ 10 times 
longer for the oscillations to be damped. Here, we neglect 
this effect, leaving its evaluation to future work. However, we 
note that it would likely result in a larger fraction of SMBHs 
remaining displaced from the nucleus, and possibly a larger 
fraction of QSOs exhibiting the correlations we propose. 
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3 TORUS MODEL 


According to the unified scheme of AGNs, galactic nuclei are 
obscured by optically and geometrically thick dusty molec¬ 
ular tori, with the amount of o bscuration depending on the 
viewing angle (I Antonuccill 1993h . Even though dust typically 
only constitutes ~ 1 per cent of the mass of the tori (while 
the rest of the mass is in gas, mostly in molecular hydrogen 
form), dust is the main component responsible for the obscu¬ 
ration and reddening of QSOs, and thus, most torus models 
focus on the dust distribution in the tori (i.e. on ‘dust tori’). 

Torus models developed to reproduce the observed spec¬ 
tral energy distributions of AGNs can be divided into two 
categor ies: smooth versus clumpy models (for an overview, 
see e.g. iHonidlioOSh . We implem ented a hydrostatic mode l 
of smooth dust tori presented bv ISchartmann et al.l (l2005h . 
This model has the advantages of simplicity, using physically 
reasonable assumptions about torus formation, and provid¬ 
ing good fits to the mid-infrared spectral energy distribu¬ 
tions of AGNs with relatively few free parameters. The dust 
torus is assumed to form from gas released through stellar 
winds and ejection of planetary nebulae in the nuclear star 
cluster. Even though the ejecta of individual stars should 
produce a cloudy structure of the torus, no instruments so 
far has been able to resolve single clo uds of the dust dis¬ 
tribution, and lSchartmann et al.l (|200^ conclude that these 
clouds should be small. For simplicity, they treat the torus as 
a continuous medium characterized by the density distribu¬ 
tion pd (see their equation 8). The effective potential created 
by the central SMBH and the angular momentum distribu¬ 
tion in the stellar core makes pa axisymmetric around the 
galactic centre. 

To cover the range o f masses of the nuc lear star 
clusters in the examples in Schartmann_et_^ ([200^ and the 
observations presented in lLeigh. Boker fc Knigga l 2012ll . we 
chose M* = where P was drawn from a uniform 

distribution in the range /3 £ [0,4]. We set the mass of the 
dust enclosed in the torus to Mdust = (5.79 x 10 ®Mq) x 
[M*/(2 X 10® Mq)] in order to reproduc e parameter values 
of the example torus given in table 1 in ISchartmann et aP 
(|2005[1 . 

The radius of the stellar core, 7?c, was chosen to be 


i?c = 2 X pc, 


( 1 ) 


where the dimensionless constants a = 0.3042 and b = 
1.2679 were obtained by fitting the half-light core radii, 
^ef f = 0.5Rci for late - and early-type galaxies in Figure 13 
in iGeorgiev fc Bokeil ll2014ll . Note that we assumed a con¬ 
stant stellar mass-to-light ratio, Mt/Lv, where Lv is the 
total luminosity of the nuclear star cluster in the V band. 

The radius Rt of the torus was chosen such that the 
sizes of the torus and of the nuclear star cluster are compa- 
rable, with Rt = 5 pc f or M* = 2 x 10 ®Mq (see table 1 in 
ISchartmann et ^l2005l l: 

Rt = 3.56 X 10“® Rc + 2.33 pc. (2) 

For all other param eters in our torus s imula tions, we 
followed the methods in ISchartmann et al.l (l2005l l. Specifi¬ 
cally, we set the outer radius of the torus Rout = 377c, the 
exponent of the angular momentum distribution in the stel¬ 
lar core 7 = —0.5, and the turbulent velocity of the clouds 
building up the torus Ut ~ cr*, where we used equation (9) in 
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Figure 3. An example for the adopted geometry and interior 
density structu re along a cross-section o f a torus surrounding an 
AGN, based on ISchartmann et al.l ll2005l V The mass of the nuclear 
star cluster containing the torus and producing its dust content 
was chosen to be M* = 10® Mq, while the mass of the central 
SMBH was set to M = 10® Mq. 


ISchartmann et ^ (|20n5h to calculate the velocity dispersion 
of the stars in the nuclear star cluster, cr«. 

The density distribution of a torus in this model is fully 
determined by two parameters: the mass of the nuclear star 
cluster, M,, and the mass of the central SMBH, M. Using 
the random values of M* and M, for each torus we calculated 
Pd{R, z) = pd{R, z)/p°d (see Eq. 8 in ISchartmann et al.ll20n5h 
in a 77out x Rout rectangular cross-section along the R — z 
plane with a linear resolution of 77out/200, set pd{R,z) = 0 
wherever Pd{R, z) < pd(Rout,0), and calculated such that 
the total mass of the torus with the resulting pa (77, z) equals 
TUdust. In our simulation, only the dust column mass density, 
Edust, was recorded as the final output, which was calculated 
by numerically integrating pd{R,z) along the line of sight 
from each SMBH position to the observer for a randomly 
oriented torus. 

For illustration. Fig. [3] shows a visualization of the ge¬ 
ometry and interior density structure (pa) along a cross- 
section of a torus with M* = 10® Mq and M = 10® Mq. 
Additionally, in Fig. [I] we show the histogram of the outer 
radii, 77out, of the 2 x 10® dust tori. This shows that the 
outer radii follow a oc l/77out distribution, ranging from 
77out,min — 10 pC tO 77out,max — 180 pC. 

4 RESULTS 

We are now ready to present the results of the Monte Carlo 
simulation of 2.5 x 10® random observations of recoiling 
QSOs. In Fig. [5] we show the histogram of dust column 
mass densities (Edust) obtained by integrating the mass den¬ 
sity of the randomly oriented dust tori along the line-of-sight 
to the recoiling SMBH (dark blue curve). For reference, we 
also show a Edust histogram for a second set of simulations 
of 2.5 X 10® QSOs with recoil velocities set to Urecoii = 0 
(fight red curve). We have found that the total number of 
unobscured QSOs (i.e. with Edust = 0) is ~ 152, 700 (~ 61 
per cent) and ~ 44, 200 (~ 18 per cent) for the Urecoii > 0 
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Figure 4. Histogram of outer radii, Rout , of the 2 x 10^ simulated 
dust tori. -Rout values follow a oc 1/Rout distribution, ranging 
from Rout ,min ~ 10 pc to Rout ,max ~ 180 pc. 



Figure 5. Histogram of dust column mass densities (Sdust) for 
2.5 X 10® random observations of QSOs, obtained by integrating 
the mass density of the randomly oriented dust tori along the 
line of sight to the recoiling, off-centre SMBHs (dark blue curve). 
For reference, we show the same histogram for a similar sample of 
QSO without any recoil (Urecoil = 0; light rod curve). The recoiling 
QSOs are much more likely to be seen completely unobscured 
(there are ~ 152700 (~ 61 per cent) and ~ 44200 (~ 18 per cent) 
QSOs with Sdust = 0 in the samples with and without recoil, 
respectively). The recoiling QSOs also show a tail of high Sdust 
values up to twice the maximum value for non-recoiling QSOs, as 
expected (see text). 


and Wrecoii = 0 samples, respectively. QSOs in the Vrecoii = 0 
sample do not have Edust values above Edust — 6 g/m^ be¬ 
cause such QSOs are never obscured by more than a half 
cross-section of thei r dust torus, and a half c ross-section of 
any dust tori in the lSchartmann et al.l (l2005l ) model have a 
maximum possible dust column density of Etorus — 6 g/m^. 
Recoiling QSOs (Urecoii > 0 km/s) can obtain up to twice 
higher values of Edust (i.e. up to Edust — 12 g/m^), which 
corresponds to the maximal obscuration of a QSO located 
in the ‘equatorial plane’, behind the dust torus (and for the 
geometrically largest tori). This configuration is rare, with 
only 409 (~ 0.2 per cent) QSOs in the UrecoU > 0 sample 
having Edust > 6 g/m^. 

Fig-[6]shows Edust versus | An| for the 2.5 x 10® recoiling 
QSOs. We have calculated the sample correlation coefficient, 
r, between Edust and logjQ(|Aii|), defined by 

^ __ (Edust logiod An|)) — (Edust)(logio(| An|)) _ 

V(Sdust) - (Sdust) V(logio(l^«l) > - (logio(|An|)}2’ 

(3) 

restricted to various ranges of |At;|niin ^ |An| ^ |An|max- 
Here (...) refers to averaging over the sample of 2.5 x 10® 
QSOs, or its various subsets. 


Fig. [3 shows r as a function of |At)|niin (left-hand pan¬ 
els) and |Au|niax (right-hand panels) for sub-samples of ob¬ 
scured QSOs with |Au| ^ |Au|min and |Aw| < |All|inax, 
respectively. As can be seen in the figure, obscured QSOs 
with |Au| < 5 km/s show no correlation between their Edust 
and |Aii| values (r <5 ~ 0, and the corresponding p-value 
for the hypothesis of no correlation is p ~ 0.95), while 
QSOs with I An I ^ 5 km/s show a significant correlation, 
with rs = 0.22 ± 0.01 and p ~ 0. The correlation coeffi¬ 
cient increases with |An|min until it reaches its maximum at 
I An|min = 45 km/s with r 45 = 0.28 ± 0.02. 

Allowing |An|niin and |Aw|niax to vary simultaneously 
does not change the result that r is the highest for the 3,824 
obscured QSOs with | Au| ^ 45 km/s. Depending on the sig¬ 
nal to noise, the measurement error on |A?;| of individual 
QSOs is a few tens of km/s (for example from spectral line 
fitting, or by using the centroids of indi vidual broad lines; 
see e.g., I.Tu et al.ll201,3l : IShen et ^l2013l ). We therefore re¬ 
strict our further investigations to the 3824 obscured QSOs 
with |Aw| ^ 45 km/s. For this subset, we have found the ap¬ 
proximate relation Edust (log,Q{| Au|}) = (1.7 ± 0.1) g/rsA x 
(logroll Au|/[1 km/s]}) “ (1.3 ± 0.2) g/m^, although of 

course this relation depends strongly on the underlying set 
of models. We also note that in addition to the random 
measurement errors on |Aw|, broad line centroids can shift 
by 0(100 km/s) between pairs of observati ons only weeks 
apart, for physica l reasons unrelated to recoil il.Tu et al.ll2013l : 
IShen et al]l2013l ). These shifts, as well as the measurement 
errors, are not expected to correlate with Edust- However, 
they introduce random errors on the Edust (log^QdAu|)) re¬ 
lation, and can therefore increase the number of QSOs re¬ 
quired to detect the underlying Edust “logiQ(|Au|) correla¬ 
tion (see discussion below). 

We next estimate the minimum number of QSOs, Amin, 
required for a significant detection of the correlation. To 
do this, for each trial N G [10,350] QSOs, we (i) gen¬ 
erated 1000 independent random samples of N QSOs, se¬ 
lected from among the 3824 obscured QSOs recoiling with 
I Au| 45 km/s, (ii) computed the p-value of the hypothesis 
of no correlation between Edust and |An| for each of these 
1000 samples, and (iii) measured the fraction of samples for 
which p < 3 X 10“® (corresponding to the rejection of the 
no-correlation hypothesis at 3cr confidence). This fraction, 
Fsct, is shown in Fig. [8] (blue crosses, right-hand panel), to¬ 
gether with the average (r) over the 1000 random samples, 
as a function of N (blue crosses, left-hand panel). 

As a test of these statistics, we have repeated the cal¬ 
culations of both (r) and F^cr for two additional sequences 
of 1000 random samples of A G [10, 350] QSOs. In the first 
set, we began with recoiling QSOs with |Au| < 5 km/s; the 
results are shown by the red empty circles in Fig. mi- In the 
second set, we computed random Edust values for stationary 
SMBHs (urecoii = 0), and associated a |Au| with each case, 
drawn randomly from the |Au| distribution for the recoiling 
(urecoii > 0) SMBHs. The results for these ‘scrambled’ data” 
are shown in Fig. |8]by the filled green circles. 

As can be seen in Fig. [S] neither the |An| < 5 km/s, 
nor the ‘scrambled data’ sets show correlations, i.e. (r) ~ 0 
and Ficr — 0 for all A. On the other hand, for obscured 
QSOs with |Au| ^ 45 km/s, we find (r) ~ 0.28 for all As, 
and Fzcr increases monotonically with A from F^a ~ 0 to 
Fsct — 1 within the interval A G [10,350]. We conclude 
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Figure 6. The line-of-sight dust column density, S^usti versus 
the (absolute value) of line-of-sight velocity offset |Au| relative 
to the host galaxy, for each of the 2.5 X 10® recoiling QSO. For 
completeness, we show the entire range of velocities we simulated, 
down to unmeasurably low values. To illustrate how far each re¬ 
coiling QSO has travelled from the galactic centre by the time 
of the mock observation, we sorted the data into three groups, 
based on their radial offset d from the galactic centre: QSOs with 
the smallest offset (d Sf 5 pc; blue crosses), QSOs with offsets 
within roughly an order of magnitude of the size of their par¬ 
ent torus (5 pc < d < 1 kpc; red circles), and the most dis¬ 
tant QSOs (d ^ 1 kpc; green full circles). The distribution of 
points shows three distinct Sdust values: Sdust = 0 (correspond¬ 
ing to the ejected QSOs with the largest offsets, viewed in an 
unobscured direct ion), Sdnat — 6 g/m^ (th e maximum possible 
obscuration in the ISchartmann et al.l ll2005h model by half of the 
torus, for QSOs whose orbits have fully decayed back to the centre 
by dynamical friction and are viewed along the equator), and at 
^dust — 12 g/m^ (the maximum possible obscuration, for QSOs 
located behind the full torus). 



Figure 7. Sample correlation coefficient r between Sdust a-ixi 
logiodAt;|), with its la uncertainty (upper panels, blue curves), 
and the corresponding p values for the hypothesis of no correlation 
(lower panels, red curves). The left-hand (right-hand) panels show 
both r and p as a function of |Au|min (|Au|niax), for sub-samples 
of obscured QSOs with |Au| ^ |Au|min (|Ati| < |Au|max). QSOs 
with I At) I > 5 km/s show a positive correlation, with r as high as 
r = 0.28 ± 0.02 for | Ai)|niin R) 45 km/s. 


that the hypothesis of no correlation between Edust and 
logic (I At)I) could be rejected with a 3cr significance 95 per 
cent of the time by observing a random sample of ~ 260 
obscured QSOs with |At)| 15 45 km/s. Note that in our sim¬ 
ulations, 3824 QSOs (~ 1.5 per cent) were obscured having 
|At)| ^ 45 km/s, thus gathering a sample of ~ 260 such 
QSOs may require observing as many as ~ 17000 recoiling 
QSOs (and a larger number of QSOs in general, accounting 



Figure 8. Left-hand panel: the expectation value (r) of the cor¬ 
relation coefficient between Sdust and logiQ(|Au|), computed in 
1000 random samples of N QSOs. The results are shown as a func¬ 
tion of N in the range N S [10, 350] for three different QSO data 
sets: obscured QSOs with | An] ^ 45 km/s (blue crosses), obscured 
QSOs with |Au| < 5 km/s (red circles), and QSOs without recoil, 
but with I At) I values assigned randomly from the |At)| distribu¬ 
tion for recoiling SMBHs (‘scrambled data’; filled green circles). 
While (r) ~ 0 tor all N in the last two data sets, the first data set 
shows a stable correlation with (r) ~ 0.28 for all N. Right-hand 
panel: the fraction F^a of the 1,000 random samples in which the 
‘no correlation’ null hypothesis can be rejected at ^ 3a signifi¬ 
cance. While Fdo- ~ 0 for the | Au| < 5 km/s and ‘scrambled’ data 
sets, it increases monotonically with N from ~ 0 to 1 for the 
[Auj 15 45 km/s data set, rising above Fso- 15 0.95 for N 15 260. 


for the fact that not all QSOs experienced a recoil in the 
past). 

To examine the impact of the sign of the line-of-sight 
velocity offset, in Fig. [9] we show Eduat versus |At)| in the 
subsample of 17157 obscured QSOs with |At)| ^ 5 km/s, 
and indicate whether they are moving away from (At) > 0, 
shown as blue crosses) or towards (At; < 0, shown as red 
full circles) the observer. First, we note that receding QSOs 
are slightly more common, with an excess of 1405 and 1162 
in the |At)| ^ 5 km/s, and |At)| ^ 45 km/s samples, respec¬ 
tively, over the approaching cases. As Fig.^shows, Edust val¬ 
ues of the receding QSOs also tend to be higher, especially 
when the velocity offset is large (|At)| > 100 km/s). This is 
because receding QSOs are more likely to be found behind, 
rather than in front of the tori. In order for a SMBH to be 
approaching with a large line-of-sight velocity while being 
located behind the torus, it would have had to reach its first 
U-turn. As mentioned above, the majority of SMBHs with 
large recoil velocities do not make their first U-turn within 
a QSO lifetime (or within our mock observation times). 

In Fig. 1101 we compare the numerically constructed 
probability density functions of various parameter distribu¬ 
tions of QSOs within the whole sample, and within the sub¬ 
sample of obscured QSOs with |Au| [5 45 km/s. Based on 
the comparison, we conclude that AGNs where the central 

SMBH has a mass around M ~ (0.5-3) x 10 ®Mq, a 

positional offset from the centre of the host galaxy between 
dofiset — 10“^ —10^ pc, and having a nuclear star cluster with 
mass M* > Mq, have the highest probability of being in 
the subset of obscured QSOs with | Au| ^ 45 km/s. The fact 
that high-mass nuclear star clusters are over-represented in 
the subset of QSOs showing Edust~|At)| correlation (see the 
upper-right panel in Fig. [13 suggests that such QSOs are 
more likely to be found in active galaxies with higher masses. 
We attribute this to the fact that in more massive galaxies. 


© 2015 RAS, MNRAS 000, HU] 








































A method to search for recoiling SMBHs in AGNs 7 


c\j 



H 



10 ^ 10 ^ 10 ^ 
lAvl [km/s] 



Figure 9. Sdust versus |Ai;| for the 17157 obscured QSOs with 
|At)| ^ 5 km/s. Receding (i.e. Au > 0) QSOs are shown by 
blue crosses, while approaching QSOs (i.e. An < 0) are shown 
by red filled circles. We also show the E^ust (logiod^'^D) ~ 
(1.7 ±0.1) g/m2 X (logio(|Ai)|/[l km/s])) - (1.3 ±0.2) g/m^ 
curve (green solid line) we obtained by fitting E^ust versus 
logio(|Au|) values for obscured QSOs with |Au| ^ 45 km/s. Re¬ 
ceding SMBHs are slightly more common among the obscured 
(^dust > 0) QSOs, their velocity offsets extend to higher values, 
and they tend to be more heavily obscured (larger Edust)) espe¬ 
cially at large velocity offsets (An > 100 km/s). See text for the 
origin of these differences. 


Figure 11. The fraction of obscured (Edust > 0) sources 
among the entire sample of 2.5 x 10^ QSOs. The trend of 
^obs(|Au|) decreasing with |Au| from .7^obs(< 10 km/s) > 0.8 to 
^obs(> 10^ km/s) < 0.4 could provide an independent test for 
a chosen combination of recoil, trajectory, and dust tori mod¬ 
els, when compared with the observed fraction of obscured (e.g. 
type II) QSOs in the future. 


observed fraction of obscured (e.g. type II) QSOs, which can 
potentially provide an independent test for a chosen combi¬ 
nation of recoil, trajectory, and dust tori models. 






Figure 10. The probability distribution functions of different pa¬ 
rameters of the QSOs within the entire sample (blue curves), and 
within the subsample of obscured QSOs with significant correla¬ 
tions (i.e. with |Ai;| ^ 45 km/s; red curves). The four different 
panels, clockwise from top left, show the SMBH masses (M), the 
masses of the nuclear star clusters (M*), recoil velocities (urecoil)? 
and the offsets of QSOs from the centres of their host galax¬ 
ies (^offset)- The difference between the red and blue curves re¬ 
veal that QSOs with M ~ (0.5 - 3) x 10® Mq, M* > IO^Mq, 

'^recoil — 100 km/s, and ^offset — 10“^ — 10^ pc are the most 
over-represented in the sample showing strong correlations be¬ 
tween Edust ^i^d |Ai;|. 


the recoiling BHs spend longer times at off-centre distances 
comparable to the size of the obscuring torus. This host-size 
dependence of the proposed correlation is another potential 
testable prediction of our toy model. 

Finally, we have studied the fraction of obscured QSOs 
within the whole sample of QSOs, J^obs(|Au|). We have 
found that this fraction decreases monotonically with |Au| 
from JTbs(< 10 km/s) > 0.8 to JTbs(> 10^ km/s) < 0.4 (see 
Fig. ini). This result can be compared in the future to the 


5 DISCUSSION 

The results in the previous section show that a positive cor¬ 
relation may be present between Edust and |Au|. We em¬ 
phasize that this conclusion is based on a specific set of 
idealized toy models - more realistic models will undoubt¬ 
edly lead to different predictions. Nevertheless, our results 
should motivate an observational search; more realistic and 
flexible models will likely be needed to interpret any correla¬ 
tion discovered in actual data. We plan to present the results 
from a search in the SDSS DRIO data base in a forthcoming 
publication. 

Here we address two natural caveats: while proxies can 
be used to estimate both Edust and |Au|, these estimates will 
have uncertainties we have so far ignored. Although we have 
no a priori reason to expect that these uncertainties will be 
correlated, they will represent noise, which will reduce the 
statistical significance of any correlation, and thus increase 
the number of QSOs required for a detection. Furthermore, 
we have so far assumed that a clean QSO sub-sample with 
Edust > 0 and |Au| ^ 45 km/s can be constructed. More 
realistically, observational errors will pollute any such selec¬ 
tion (e.g. QSOs with |Au| ^ 45 km/s can be missed, and 
QSOs with |At)| < 45 km/s can be mistakenly included). 

The line-of-sight velocity offset (|At;|) can be estimated 
using the centroids or peaks of broad lines, relative to those 
of the narrow lines. The level of obscuration (Edust) can 
be estimated using QSO colours: a larger obs curing column 
should result in a redder continuum (see e.g. iLedoux et al.l 
ll2015l) and references therein). An alternative proxy for the 
obscuring column could be the equivalent width of the broad 
lines themselves. Although the geometry of the broad line 
region is not fully understood, type II quasars have weak 
or no br oad emission lin es. In the standard unified model 
of AGN llAntonuccill 19931) . this is attributed to the obscura¬ 
tion of the high-velocity broad-line region by the dusty torus 
(note that the continuum is still visible and can therefore not 
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be spatially colocated with the broad-line emitting region). 
This implies that the equivalent width of these broad lines 
must vary with the level of the obscuration. 

Here we estimate the impact of having random errors on 
both |Aii| and Edust. As mentioned above, we expect |Au|, 
when determined using line centroids, to suffer physical vari¬ 
ations of order aAv ~ 100 km/s. We therefore repeated our 
analysis above, except we added a random additional An 
component to the velocity offset of each simulated QSO, 
drawn from a Gaussian distribution with a standard devia¬ 
tion of aAv For any given ctad, we then selected the subset 
of the 2.5 X 10® QSOs with Edust > 0 and |An| ^ 45 km/s, 
as before. Finally, we measured the expectation value of the 
correlation coefficient (r) between Edust > 0 and |An| in 
1000 randomly drawn sets of N QSOs from this subset, and 
determined the minimum number Nmin of QSOs required to 
infer r > 0 with 3a significance 95 per cent of the time. 

The result of this random-velocity-error exercise is 
shown in Fig. [T5] for both Edust-logj^d Aii|) and Edust-|Aii| 
correlations. The right-hand panel in this figure shows Amin 
as a function of the random error a^v As expected, the 
number of QSOs required to detect a correlation increase 
monotonically with (tau, rising to Amin = 3560 for the 
Edust-|An| correlation and for aAv = 100 km/s; this is more 
than an order of magnitude increase over the Amin = 260 
in the idealized case for the Edust-logjQ(|At)|) correlation 
without any velocity-offset error. The left-hand panel in 
Fig. [12] shows the corresponding values of the correlation 
coefficients; (r) decreases monotonically with aAv, reaching 
(r) = 0.091 ± 0.016 for the Edust-|An| correlation and for 
aAv = 100 km/s. This is a decrease from the (r) = 0.28±0.02 
in the idealized case for the Edust-logjQ(|At)|) correlation by 
a factor of three. Additionally, we have calculated (r) for 
partially obscured QSOs only (i.e. QSOs with 0 < Edust < 
2.3 g/nF), again assuming aAv ~ 100 km/s, and we have 
found it to be (r) = 0.026 ± 0.004, measured between |At;| 
and Edust- 

Similarly to the errors in |At;|, we have studied the 
impact of errors in Edust. Here, there are several poten¬ 
tial concerns. First, in the presence of uncertainties, a 
sharp selection Edust > 0 is unrealistic. Secondly, for large 
obscuring columns, with Edust comparable to the maxi- 
mu m possible column densit y of a half-torus according to 
the ISchartmann et ^ (l2005l 'l model (see Fig. O , Etorus — 
6 g/m^ (i.e. for configurations corresponding to true type 
H QSOs), the broad lines may be undetectable, or may 
have too low S/N for a reliable centroid measurement □ 
We here repeated our fiducial analysis, except that we re¬ 
placed the selection criterion Edust > 0 by the range Emiu < 
Edust ^ Emax- A^e chose Emin — 0.6 g/m uS 0.1 X Etorus and 
Emax = 2.3 g/m^ « 0.4x Etorus- The number of QSOs within 
this range (and also satisfying Av > 45 km/s as before) is 
A = 1436, or ~ 0-6 per cent of the full sample of 2.5 x 10® 
QSOs. With this cut, we find that the correlations remain 
significant: (r) = 0.098 ± 0.026 for the Edust-logtodAu|) 
correlation, with a p-value of p = 2 x 10“^. 

We also evaluated the impact of random Gaussian errors 


^ Note that the broad lines and the optical continuum are 
thought to arise from different spatial locations, and can suffer 
different levels of obscuration; we ignore this complication here. 



Figure 12. The figure shows the impact of random errors in the 
velocity offset; the error is assumed to be Gaussian with a stan¬ 
dard deviation of ctav The right-hand panel shows the increase 
in the minimum number of QSOs (Amin) required to detect cor¬ 
relations between Edust a-nd log]^g(|Au|) (blue solid curves), or 
between Edust and |Au| (red dashed curves) in a subset of QSOs 
with Edust > 0 and |Au| ^ 45 km/s, as uav increases. Curves in 
the left-hand panel, with ilcr errors, show the corresponding de¬ 
crease in the correlation coefficient (r) (see text for more details). 
Note that measurement errors below (tau ^ 45 km/s results with 
higher (r) and lower Amin values for the EdustAog]^Q(|Au|) cor¬ 
relation, while the Edust^|Au| correlation appears to be stronger 
when (TAu ^ 45 km/s. 


in Edust > 0 with a standard deviation of as. The results are 
similar to those shown in Fig. 1121 we find that (r) decreases, 
and the minimum number of QSOs required to detect the 
correlation between Edust and logj^gd At;|) (Amin) increases 
with as- For the case of as ~ 0.6 g/nF = 0.1 x Etorus, we 
find (r) = 0.27 ± 0.05, with a p-value of p « 10“®®, and 
Amin = 290. Thus, random Gaussian errors in Edust > 0 
have a very small effect on the Edust“logjQ(|Au|) correlation. 

Based on the above, we conclude that statistical errors 
on estimating the velocity offsets and obscuring column den¬ 
sities can diminish the correlations predicted in the idealized 
toy model. This is mostly due to the impact of these errors 
on the purity of selection of a sub-sample, intended to in¬ 
clude only the most correlated QSOs. In the examples above, 
we find that this increases the number of QSOs required for 
a detection by an order of magnitude, but there is no in¬ 
dication that these errors will prevent a detection of the 
correlations in a large sample (tens of thousands) of QSOs. 

Finally, as the distribution of SMBH masses in QSOs 
observable in different survey projects may vary signifi- 
cant ly, we have considered how usi ng higher SMBH ma sses 
(see iGreene fc h 3|2007I. l2009l and IShen fc Kellvll2012ll af¬ 
fect our results. In these tests, we have also taken into ac¬ 
count the possibility of a correlation between the SMBH 
mass and the mass of the stellar core llFerrarese et al.l[2006l : 
iGraham et al.l[201lj : iKormendv fc Holl2013lf . We have found 
that the increased SMBH masses lead fewer SMBHs escap¬ 
ing the galactic nucleus (an effect that should increase the 
number of SMBHs in the population showing Edust-|An| 
correlation), however they also lead to significantly smaller 
torus sizes, resulting with a net decrease in the correlation 
coefficient of a factor of ~ 2. This demonstrates that al¬ 
though the strength of the Edust“| An| is sensitive to the un¬ 
derlying assumptions of the applied models, the Edust-|An| 
correlation itself is a robust result. 
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6 CONCLUSIONS 

A large fraction of luminous QSOs are believed to be trig¬ 
gered by merger events, and if QSO activity often follows 
the coalescence of SMBHs in the merging nuclei, then a large 
fraction of all QSOs may have experienced recent gravita¬ 
tional recoil. While a handful of QSOs with large spatial 
or velocity offse ts have been identified as recoil candidates 
llKomossall2012l 'l. the ubiquity of this phenomenon remains 
poorly known. 

Here, we proposed a statistical technique to search for 
a population of recoiling SMBHs with smaller velocity off¬ 
sets among luminous QSOs, based on a positive correlation 
between these velocity offsets and the column densities of 
obscuring dust tori. The correlations are introduced by the 
damped oscillating motions of SMBHs with typical recoil 
speeds, which are comparable to the escape velocities from 
galactic nuclei. These SMBHs spend a significant time off- 
nucleus, and at distances comparable to the size of obscuring 
tori in located in the nuclei. 

We have demonstrated, using simple models of gravi¬ 
tational recoil, SMBH trajectories, and the geometry of ob¬ 
scuring dust tori, that observing a random sample of a few 
thousand partially obscured QSOs, with line-of-sight veloc¬ 
ity offsets |Au| ^ 45 km/s, could allow a significant detec¬ 
tion of the correlation between the line of sight dust column 
density Edust and |Au| (or logio(|Au|)). 

The existence of this correlation is found to be robust 
within the simple models we choose, and detectable even in 
the face of random errors on Edust and |Au|. However, we 
expect that the tightness of the correlation, the underlying 
Edust — |Au| relation, and the properties of the QSO subset 
showing correlations can depend strongly on the model as¬ 
sumptions. Nevertheless, our results should motivate search¬ 
ing for correlations in real data; a positive detection would 
allow testing various combinations of recoil, BH trajectory, 
and dust tori models. A further test could be possible by 
comparing the predicted fraction of obscured QSOs to the 
observed fraction of type H QSOs. We propose to carry out 
these tests in the future using catalogues of observational 
data for QSOs, such as the SDSS-DRIO QSO catalogue. 


ACKNOWLEDGEMENTS 

The authors would like to thank Bence Becsy, Gergely 
Dalya, and Akos Szolgyen for useful discussions and for valu¬ 
able comments on the manuscript. Peter Raffai is grateful for 
the support of the Hungarian Academy of Sciences through 
the ’Bolyai Janos’ Research Scholarship programme. We also 
acknowledge financial support from the NASA Astrophysics 
Theory Program under grant no. NNX11AE05G (to ZH) 
and from OTKA under grant no. 101666 (to ZF). 


Begelman M. C., Blandford R. D., & Rees M. J., 1980, 
Nature, 287, 307 

Blecha L. et ah, 2011, MNRAS, 412, 2154 
Blecha L. et ah, 2015, preprint (arXiv: 1508.01524) 
Bogdanovic T., Reynolds C. S., & Miller M. C., 2007, ApJ, 
661, L147 

Bonning E. W., Shields G. A., & Salviander S., 2007, ApJ, 
666, L13 

Gomerford J. M. et ah, 2009, ApJ, 698, 956 
Dubois Y. et ah, 2014, MNRAS, 440, 2333 
Ferrarese L. et ah, 2006, ApJ, 644, L21 
Georgiev I. Y., & Boker T., 2014, MNRAS, 441, 3570 
Gerosa D. et ah, 2015, MNRAS, 451, 3941 
Graham A. W. et ah, 2011, MNRAS, 412, 2211 

Greene J. F., & Ho L. G., 2006, ApJ, 641, L21 

Greene J. F., & Ho L. G., 2007, ApJ, 667, 131 

Greene J. F., & Ho L. G., 2009, ApJ, 704, 1743 

Gualandris A., & Merritt D., 2008, ApJ, 678, 780 
Guedes J. et ah, 2009, ApJ, 702, 890 
Guedes J. et ah, 2011, ApJ, 729, 125 
Haehnelt M. G., 1994, MNRAS, 269, 199 
Healy J., Lousto C. O., & Zlochower Y., 2014, Phys. Rev. 
D, 90, 104004 

Honig S. F., 2008, Dr. rer. nat. dissertation, Rheinische 
Friedrich-Wilhelms-Universitat 
Ivanov P. B., Papaloizou J. C. B., & Polnarev A. G., 1999, 
MNRAS, 307, 79 
Ju W. et ah, 2013, ApJ, 777, 44 
Komossa S., 2012, Adv. Astron., 2012, 364973 
Komossa S., & Merritt, D., 2008, ApJ, 689, L89 
Kormendy J., & Ho L. C., 2013, ARA&A, 51, 511 
Ledoux G. et ah, 2015, A&A, 580, A8 
Leigh N., Boker T., & Knigge G., 2012, MNRAS, 424, 2130 
Loeb A., 2007, Phys. Rev. Lett., 99, 041103 
Lousto C. O. et ah, 2012, Phys. Rev. D, 85, 084015 
Madau P., & Quataert F., 2004, ApJ, 606, L17 
Martini P., 2004, in Ho L. G., ed., Coevolution of Black 
Holes and Galaxies, Cambridge Univ. Press, Cambridge, 
p. 170 

Paris 1. et ah, 2014, A&A, 563, A54 
Reynolds C. S., 2013, Class. Quantum Gravity, 30, 244004 
Schartmann M. et ah, 2005, A&A, 437, 861 
Sesana A., Haardt F., Madau P., & Volonteri M. 2005, ApJ, 
623, 23 

Shen Y., & Kelly B. C., 2012, ApJ, 746, 169 
Shen Y., Liu X., Loeb A. & Tremaine S., 2013, ApJ, 775, 
49 

Tanaka T. & Haiman Z., 2009, ApJ, 696, 1798 

Tremaine S. et ah, 2002, ApJ, 574, 740 

Volonteri M., Haardt F., & Madau P., 2003, ApJ, 582, 559 


REFERENCES 

Aller M. C., & Richstone D., 2002, AJ, 124, 3035 
Antonucci R., 1993, Ann. Rev. A. & A., 31 473 
Baker J. G. et ah, 2008, ApJ, 682, L29 
Barth A. J., Greene J. E., & Ho L. C., 2005, ApJ, 619, 
L151 


© 2015 RAS, MNRAS 000, [TH9] 




